function out = genclutter(beta_FA, range)

V   = (range(1,2) - range(1,1))*(range(2,2) - range(2,1));
m_k = poissrnd(beta_FA*V);
x   = rand(1,m_k)*(range(1,2) - range(1,1)) + range(1,1);
y   = rand(1,m_k)*(range(2,2) - range(2,1)) + range(2,1);

out = [x; y];